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Abstract 

Scalar cosmological perturbations with nearly flat power spectrum may originate 
from perturbations of the phase of a scalar field conformally coupled to gravity and 
rolling down negative quartic potential. We consider a version of this scenario whose 
specific property is a long intermediate stage between the end of conformal rolling and 
horizon exit of the phase perturbations. Such a stage is natural, e.g., in cosmologies 
with ekpyrosis or genesis. Its existence results in small negative scalar tilt, statistical 
anisotropy of all even multipoles starting from quardupole of general structure (in 
contrast to the usually discussed single quadrupole of special type) and non-Gaussianity 
of a peculiar form. 

1 Introduction and summary 

By far the most developed hypothesis on the origin of the cosmological perturbations is 
the slow roll inflation [1]. The inflationary mechanism [2] generates almost Gaussian scalar 
perturbations whose power spectrum is almost flat due to the slow evolution of relevant 
parameters (the Hubble parameter and time derivative of the inflaton field). Similar situation 
occurs in the inflationary scenario with the curvaton mechanism [3]; in either case, the 
approximate flatness of the spectrum is a direct consequence of the approximate de Sitter 
symmetry of the inflating background. 
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In quest for an alternative symmetry behind the flat scalar spectrum one naturally turns 
to conformal invariance [4, 5]. Conformal symmetry implies scale invariance, which in the 
end may be responsible for the scale-invariant scalar spectrum. An assumption of conformal 
invariance at the time the primordial perturbations are generated is in line with the viewpoint 
that the underlying theory of Nature may have conformal phase, and that the Universe may 
have started off from, or passed through that phase. 

At the present, exploratory stage it makes sense to consider this possibility in the context 
of toy models. One such model is proposed in Ref. [4]. Besides conventional Einstein gravity 
and some matter that dominates the cosmological evolution, its main ingredient is a complex 
scalar field (j) conformally coupled to gravity. Conformal invariance implies that the scalar 
potential is quartic, while the dynamics is non-trivial if its sign is negative, 

V{ct>) = -/i^r , (1) 

where his a. small parameter. One assumes that the background space-time is homogeneous, 
isotropic and spatially flat, 

ds'' ^ a\ri){drf - d^^) . (2) 



Then in terms of the field 



X(r/,x) = a{ri)(l){ri,x) 



the dynamics is the same as in fiat space-time. One further assumes that the classical 
background field Xc is homogeneous. As it rolls down its potential V{x) — — /i^|x|^, it 
approaches the late time attractor 

^'<^'>'> = ^i^;^ • 

where rji^^ is an arbitrary real parameter ( "end of roll" ; the reason for the superscript (0) in 
notation will become clear later), and we take Xc real without loss of generality. 

The point of Ref. [4] is that the behavior of the phase^ 9 = \pl Arg in the background 
(3) is very similar to what happens at infiation to the fiuctuations of a massless scalar field 
minimally coupled to gravity (e.g., infiaton itself). The phase perturbations bQ start off as 
vacuum fiuctuations and eventually freeze out. To the leading order in /i, the resulting phase 
perturbations are Gaussian and have fiat power spectrum 

The latter property is a consequence of conformal invariance and [/(l)-symmetry <h e*"0 
inherent in the model. The phase perturbations are the source of the adiabatic perturbations 



^The normalization here is chosen for future convenience. 
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in this scenario, which proceeds as follows. At large field values, the potential V(|0|) is 
assumed to be difi^erent from (1) and to have a minimum at |0| = /o; we assume that 
/o <^ MpL (see also the discussion in Section 2.1), so that the contribution of the field to 
the effective Planck mass is always negligible. At |0| ~ /o, conformal symmetry is broken, 
the radial field |0| interacts with other fields, and its oscillations about the minimum get 
damped quickly enough. To be on the safe side, we assume that the field is a spectator 
at this and earlier stages, i.e., its energy density is small compared to the energy density 
Ptot of matter that dominates the cosmological evolution. This is the case provided that 

Ip^I ~ /^Vo' « Ptot = ^Ml^H' . (5) 

Then the decay products of the field |0| do not affect the evolution of the Universe and, 
furthermore, the perturbations of that exist before the end of rolling and disappear after 
101 gets relaxed to the minimum of do not produce substantial density perturbations 

in the Universe. 

Once the radial field |0| settles down to /o, what remains are the perturbations of the 
phase, which at this point are isocurvature perturbations. They get reprocessed into adia- 
batic perturbations at much later epoch by one or another mechanism. As an example, the 
phase 6 may be pseudo-Nambu-Goldstone field, and may serve as curvation [3, 6]. Alter- 
natively, perturbations 59 may be converted into adiabatic perturbations by the modulated 
decay mechanism [7, 8]. In either case, the adiabatic perturbations inherit the correlation 
properties from the phase perturbations (with possible additional non-Gaussianity generated 
at the conversion epoch), while the amplitude of the adiabatic perturbations is, generally 
speaking, smaller than that of the phase perturbations. In view of the latter property, we 
treat our only parameter, the coupling constant /i, as free (but small). 

The scenario cannot work at the conventional hot cosmological epoch, for the following 
reason. The vacuum state of the phase perturbations 56 is weU defined at early times 
provided that these perturbations evolve in the WKB regime, which implies 

^(^i°^ ~ "n) ^ ^ 1 early times , (6) 

where k is conformal momentum. On the other hand, the property that these perturbations 
are frozen out at late times holds if 

/j(,y(o) ^ late times . (7) 

So, the scenario requires that both of these inequalities are satisfied at conformal rolling 
stage. This can only happen if the duration of that stage in conformal time is greater than 
k"^. For conformal momenta of cosmological significance this means that conformal rolling 
lasts longer (in conformal time) than the entire hot stage until the present epoch. Thus, the 
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mechanism can only work at some pre-hot epoch at which the horizon problem is solved, at 
least formally. This is similar to most other mechanisms of the generation of cosmological 
perturbations (see, however, Ref. [9]). 

At the conformal rolling stage, the dynamics of the phase perturbations S9 is governed 
solely by their interaction with the background field (3) (as well as with the radial pertur- 
bations S\x\, see below); the evolution of the scale factor a{r]) is irrelevant. After the end 
of conformal rolling, the situation is reversed. Once the radial field |0| has relaxed to the 
minimum of the scalar potential, the phase is a massless scalar field minimally coupled 
to gravity (this is true for any Nambu-Goldstonc field [10]). Since we arc talking about a 
yet unknown pre-hot epoch, it is legitimate to ask what happens to the perturbations of the 
phase right after the end of conformal rolling. Barring fine tuning, there are two possibilities 
for the perturbations 69: 

(i) they are already superhorizon in the conventional sense at that time, or 

(ii) they are still subhorizon. 

The version (i) of the scenario has been considered in Refs. [11, 12]; in that case, the phase 
perturbations do not evolve after the end of the conformal rolling stage, and the properties 
of the adiabatic perturbations are determined entirely by the dynamics at conformal rolling 
(modulo possible non-Gaussianity generated at the conversion epoch; the latter is not specific 
to the conformal rolling scenario). To subleading orders in h, this dynamics is fairly non- 
trivial, and the resulting effects include certain types of statistical anisotropy [11] and non- 
Gaussianity [12]. 

In this paper we consider the second possibility, i.e., assume that there is a long enough 
period of time after the end of conformal rolling, at which the phase perturbations remain 
subhorizon in the conventional sense. Their behavior between the end of conformal rolling 
and horizon exit depends strongly on the evolution of the scale factor at this intermediate 
stage. In order that the fiat power spectrum (4) be not grossly modified at this epoch, 
the scale factor should evolve in such a way that the dynamics of 66 is effectively nearly 
Minkowskian. Although this requirement sounds prohibitively restrictive, there are at least 
two cosmological scenarios in which it is obeyed. One is the bouncing Universe, with matter 
at the contracting stage having super-stiff equation of state, p. It is worth noting in this 
regard that stiff equation of state is preferred at the contracting stage for other reasons [13, 14] 
and is inherent, e.g., in a scalar field theory with negative exponential potential, like in the 
ekpyrotic model [15]. It is known [16] that in models with super-stiff matter at contracting 
stage, the resulting power spectrum of scalar perturbations is almost the same as that of 
massless scalar field in Minkowski space, V{k) oc k"^. This implies that the dynamics of 
the scalar field perturbations is almost Minkowskian in these models. We discuss this point 
further in Appendix A. In tractable bouncing models like those of Refs. [17, 18, 19], our phase 
perturbations exit the horizon at the contracting stage, pass through the bounce unaffected 
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(cf. Ref. [20]), remain superhorizon early at the hot expansion epoch and get reprocessed 
into adiabatic perturbations, as discussed above. 

Similar situation occurs in another scenario suitable for our purposes, namely, "genesis" 
of Ref. [5] (see also Ref. [17]). According to this scenario, the Universe is initially spatially 
flat and nearly static, stays in this nearly Minkowskian state for long time, then its expan- 
sion quickly speeds up and eventually the conventional hot epoch begins. If our conformal 
rolling stage ends up well before the start of rapid expansion, the evolution of the phase 
perturbations is again nearly Minkowskian up until the horizon exit. 

In both scenarios the relevant range of momenta is wide, provided that /q is small enough 
(but not unrcalistically small). We discuss this point in Section 2.1. So, it is legitimate to 
approximate the evolution of the phase perturbations as Minkowskian in the time interval^ 
Tj^, — e < Tj < Tji, where rji is some time after the horizon exit, and (r^* — e) is the time when 
the radial field relaxes to the minimum of T^(|0|) and the conformal rolling stage ends. We 
set e = in what follows to simplify notations; keeping e 7^ would not change our results 
(recall that the phase perturbations are frozen out well before rj — 77*). The field (5^(x, 77*), 
determined by the dynamics at the conformal rolling stage, serves as the initial condition for 
further Minkowskian evolution from 77* to rji. Barring fine tuning, the case of interest for us 
is^ 

Kvi -??*)> 1 • 

Our purpose is to study the properties of the phase perturbations at 77 = 771, as these 
properties are inherited by the adiabatic perturbations. 

To the leading order in h, we find nothing new: the phase perturbations at = r/i are 
Gaussian and have flat power spectrum. Subleading orders in h are more interesting. A 
simple way to understand what is going on is to notice that the end-of-roU time 77*, instead 
of being a constant parameter, is actually a Gaussian random field [4], 77* (x) = 77**^^ +577* (x) 
with 577* (X h. This is due to the fact that not only the phase 9 but also the radial field 
Ixl acquire perturbations at the conformal rolling stage; after freeze out, perturbations S\x\ 
can be interpreted as perturbations 577*(x). The effect of the perturbations 5rj^ on the phase 
perturbations S6 is twofold. First, the perturbations ^77^, modify the dynamics of 66 at the 
conformal rolling stage. This property is common to both cases (i) and (ii), and we make use 
of the results of Ref. [11]. The new point is that the resulting field 56'(x, 7/*(x)) serves as the 
initial condition for the Minkowskian evolution. Second, this initial condition is now imposed 
at the non-trivial hypersurface 77 = 77* (x). This is illustrated in Fig. 1. The net result is 
that the perturbation S9{'x) at the time 771 is a combination of two Gaussian random fields 
originating from vacuum fiuctuations of the phase 9 and radial field \x\, respectively (better 

^For the reason that will become clear shortly, we drop here the superscript (0) in the notation of r?*. 
^In the opposite case, the phase perturbations do not evolve between r/* and 771 , and we are back to the 
version (i) above. 
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Figure 1: Due to the perturbations of the radial field, the evolution of phase perturbations 
proceeds in inhomogeneous background. Perturbations 66 oscillate in time at early stage 
(region I), freeze out at time t] = rjxi^) and temporarily stay constant (region II) until 
the end of conformal rolling that occurs at 77 = 77* (x). Then they evolve again, now in 
nearly Minkowskian regime (region III), until the horizon exit time rji. Later on (region IV), 
perturbations 66 are superhorizon and stay constant. 



to say, from vacuum fluctuations of imaginary and real parts of with our convention of 
real background Xc)- This leads to several potentially observable effects. 

At the level of the two-point correlation function of the phase perturbation 56'(x, rji), and 
hence of the adiabatic perturbation (, we have found two effects. The flrst one is negative 
scalar tilt 

n., - 1 : 



47r2 



We note in passing that this is not a particularly strong result, as small scalar tilt in our 
scenario may also originate from weak violation of conformal invariance at the conformal 
rolling stage [21] and/or not exactly Minkowskian evolution of 66 at the intermediate stage, 
cf. Appendix A. The second effect is the statistical anisotropy: the power spectrum has the 
form 

(9) 



V^ik)=v['\k) 1 + Q(k) 
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where "P^^^ is independent of the direction of momentum (nearly flat spectrum with small 
tilt), k = k/A; is the unit vector along the momentum and Q{k) is itself a random field, 
which depends on the direction of k only. Unlike the statistical anisotropy discussed in the 
inflationary context [22, 23, 24, 25], and also in the version (i) of the conformal roUing sce- 
nario [11], the function Q{k) contains all even angular harmonics, starting from quadrupole. 
We give here the expression for (5(k) which accounts for the quadrupole component only 
(see Section 4 for the results vahd for all multipoles) 



Q{k) ^Q-Wij [kikj--6ijj , (10) 

where Wij is a general symmetric traceless tensor normalized to unity, WijWij = 1, and the 
variance of the quadrupole component (in the sense of an ensemble of universes) is 

' 327r2 ^ ' 

Of course, the precise values of the multipoles of (5(k) in our patch of the Universe are 
undetermined because of the cosmic variance. 

Due to the interaction with the perturbations 57]*, the resulting phase perturbations 
(5^(x, r^i) and their descendant perturbations are non-Gaussian (we leave aside here the 
non-Gaussianity that may be generated at the epoch of conversion of the phase perturbations 
into adiabatic ones; our scenario is not special in this respect). Their three-point correlation 
function vanishes identically due to the discrete symmetry 9^—9 (cf. Ref. [11]), while the 
four-point correlation function has a peculiar form 

(C(k)c(k)C(k')C(k')> - 4U3 4^-3 '^^'^ + ^^^^^ + • V + 

+ (k o k') + (k ^ k') . (12) 

The leading term in (12) (unity in square brackets) is the Gaussian part, while the non- 
Gaussianity is encoded in -F/vg = 0{h'^). Note that the structure of the non-Gaussian part 
is fairly similar to that of the disconnected four-point function. Note also that Fng depends 
on the angle between k and k' only. For reasons we discuss in Section 5, the notion of 
non-Gaussianity is appropriate if the angle between k' and k is small, i.e., |k — k'| <^ 1. In 
this regime, the leading behaviour of F^g is 

^ const 

Fng = — log ^ , 

TT^ |k-k'| 

where constant in the argument of logarithm cannot be reliably calculated because of the 
cosmic variance. The logarithmic behavior does not hold for arbitrarily small |k — k'|: the 
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function Fng(^ — k') flattens out most likely at |k — k'| ~ [^(^i — V*)]~^^'^j ^^'^ certainly at 
|k — k'l ~ [k{rii — r]^)]~^. So, the parameter (r/i — is detectable in principle (but, probably, 
not in practice). 

It is tempting to speculate that the negative scalar tilt Ug — 1 ~ —0.04, favoured by 
the data [26], has its origin in the dynamics we discuss in this paper. If so, our only free 
parameter h is determined from (8), h'^ ~ 0.5, while the small amplitude of the adiabatic 

perturbations is to be attributed to the mechanism that reprocesses the phase perturbations 
into adiabatic ones. In that case the statistical anisotropy is roughly of order 1, which is 
probably inconsistent with the data. On the other hand, if one attributes the small observed 
amplitude of primordial scalar perturbations, — 5 ■ 10~^ [27], entirely to the smallness 

of h, i.e., identifies Vse with 'P^, then h"^ ~ 10^^, and the statistical anisotropy is at the level 
Q ~ 10~^, while the non-Gaussianity is probably unobservable. This gives an idea of the 
range of predictions of our model. 

This paper is organized as follows. We begin in Section 2.1 with discussing the range 
of momenta of modes under study. To make the presentation self-contained, we review in 
Sections 2.2 - 2.4 the properties of the radial and phase perturbations at the conformal rolling 
stage [11]. Phase perturbations at the end of intermediate, (almost) Minkowskian stage are 
studied in Section 3 to the first non-trivial order in h. This is sufficient for evaluating the 
statistical anisotropy and non-Gaussianity in Sections 4 and 5, respectively. The calculation 
of the tilt (8) requires the analysis of oider-h^ corrections, so we postpone it to Section 6. We 
discuss in Appendix A the properties of a massless scalar field in the contracting Universe 
filled with super-stiff matter. We present in Appendices B - D technical details of the 
calculations performed in Sections 3 and 4. 



2 Conformal rolling 

Let us review the main properties of our scalar field at the stage when it rolls down its 
potential. At this stage, the theory is described by the action 

S — Sg+m + , 

where Sq+m is the action for gravity and some matter that dominates the evolution of the 
Universe, and 

g^'^d,<f>*d,<P+^<P*<f>-V{<P) 

is the action for the scalar field we are going to discuss. Here the scalar potential is given 
by (1). We assume that the field is a spectator which does not affect the cosmological 
evolution; for this reason, mixing between this field and gravitational degrees of freedom is 
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negligible. The background metric is assumed to be given by (2). One introduces the field 
X = Ci4>, and obtains its action in conformal coordinates in the Minkowskian form, 



S[x]= J d'xdrj [v^''d,x*d.x + h'\x\'] ■ 



The homogeneous background solution Xciv) to the field equation is given by (3). Recall 
that we have chosen Xc real without loss of generality. 

2.1 Momentum scales 

Before discussing field perturbations in detail, let us consider momentum scales for which 
our scenario, outlined in Section 1, is valid. According to this scenario, conformal rolling 
stage ends up when the radial field |0| becomes of order /q. This occurs at time 77/ such that 

1 

/o. 



a{r]f )h{r]f - 77!°^) 
Hence, the shortest waves obeying (7) have present momenta 

~ IiJq , 

Oo ah do 

where gq and ah are the present value of the scale factor and its value at the beginning 
of the hot stage, respectively. On the other hand, we assume that the relevant modes are 
subhorizon right after rjf, 

^ >H{vf)- (13) 



a{r]f) 

We recall our requirement (5) and find that the longest waves obeying (13) satisfy 

— > ^ ^ ML . "(^/) . ^ (14) 
Oo Oo MpL ah Oo 

We see that the relevant range of momenta is 

ikZpL 

It is wide enough, provided that the energy scale /o is sufficiently low. As an example, for 
A;max/A;min ~ (lOkpc) " V(lOGpc) we need /o < IQ-^Mpl. 

If our mechanism is supposed to work at contracting stage in the bouncing Universe 
scanario with the hot epoch starting immediately after bounce, the inequality (14) implies 
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much stronger bound on /q. Indeed, a{rif )/ah > 1 in this scenario, while ah/ do ^ T^/Th > 
To/MpL. We require that /cmin/ao is lower than the present Hubble scale Hq and obtain 

/TT \ 1/2 

• /o < MpL ( ~ lO-^^MpL . (15) 

Even for /i ~ 10~^ this imphes /o < lO^GeV. Interestingly, fully consistent with this scenario 
is the scale /o ~ TeV. 

On the contrary, a{r]f)/ah can be large in the "genesis" scenario [5]. Therefore, no bound 
similar to (15) can be established in that case. 



2.2 Radial perturbations 

Let us consider perturbations about this background. To the leading order in h, perturba- 
tions Sxi — V^S{Rex) and 6x2 — V^lmx decouple from each other. At early times, when 
the relation (6) is satisfied, the fields Sxi and 6x2 are free and Minkowskian. The vacuum is 
well defined and we assume, as usual, that this vacuum is the initial state. The normalization 
factor is chosen in such a way that the real fields 5xi and 6x2 are canonically normalized 
at early times. At late times, when the opposite inequality (7) holds, the perturbations no 
longer oscillate. 

We begin with the radial perturbations ^xi. They obey the linearized field equation, in 
momentum representation, 

{Sxi)"+P' Sxi - 6h\lSxi ^ {Sxi)" + P' Sxi - -To^—JXi = , (16) 

[r]* - V) 

where prime denotes the derivative with respect to conformal time. We denote the conformal 
momentum of the radial perturbation by p and reserve the notation k for the conformal 
momentum of the phase perturbation. The properly normalized solution to Eq. (16) is 




(0) _ 

SXi = i^\l ^-^Hi% [p(r;l°) - V)] . Bp + h.c. , 

where Bp, B^ are annihilation and creation operators obeying the standard commutational 
relation [Bp,B^,] = 6{p — p'), H^^^^ is the Hankel function, and here and in what follows 
in this Section we omit irrelevant phase factors. At late times the solution approaches the 
asymptotics 

The interpretation of the behaviour 5xi oc (77**^^ — rj)"^ is that the end-of-roU parameter 
77* becomes a random field. Indeed, with perturbations included, the radial field Rex = 
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+ 5xi/\/2 can be written at late times as follows, 



where 

and the linearization in ^77* is understood. The field ^77* is constant in time and is given by 

Note that this field has red power spectrum, 

^ ■ ^^^^ 

Clearly, the overall spatially homogeneous shift of the end-of-roll time is irrelevant, as it can 
be absorbed into redefinition of the bare parameter 7]i^\ What is important is the gradient 
of 77* (x), as well as higher derivatives. It is convenient to introduce the notation 

It reflects the fact that to the first order in the gradient expansion of ?7*(x) (i.e., neglecting 
the second derivatives of ?7*(x)) and to the linear order in v, the hypersurfaces of constant 
Rex, i.e., hypersurfaces 77* (x) — rj — const, are boosted with the velocity v with respect to 
the cosmic frame: these are hypersurfaces 77 + vx = const. The random field v(x) has flat 
power spectrum, while higher derivatives of 77* (x) have blue spectra. 



2.3 Phase perturbations: order v 

Let us now turn to the perturbations 5x2 of the imaginary part, and account for their 
interaction with radial perturbations. As we will see in what follows, relevant perturbations 
577* have wavelengths much longer than the wavelengths of the phase perturbations, 

where, as before, p and k are conformal momenta of radial and phase perturbations, respec- 
tively. Because of this separation of scales, it is legitimate to use the expression (17), valid in 
the late-time regime ^(77* — 77) ^ 1, when considering the dynamics of 5x2, and treat the field 
(17) as the background. It is worth noting, however, that the expression (17) is valid to the 
linear order in ^7;* only; furthermore, there are corrections to (17) of order didjr]^.{x) / {rj^ —rj) . 
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Therefore, the resuhs of this Section are vahd to order v (or, equivalently, to the subleading 
order in h). We present the expressions vahd to order v"^ in Section 2.4. 
With this quahfication, the hnearized field equation for 8x2 reads 

{5x2)" - OA Sx2 - ^h^Rexf ■ 6x2 = (6x2)" - dA Sx2 - . , f_ J X2 = . (19) 

At early times, when fc(r/* — r/) S> 1, we get back to the Minkowskian massless equation, and 
the solutions are spatial Fourier modes that oscillate in time. Hence, the solution to Eq. (19) 
has the following form, 

= / (27r)t^v^ ('^Xrnk,x,r,)i. + /..c.) , 

where 6x2 ^(k, x, ?]) tends to q^^^^'-^'^ 7^ _qo A^, A\. is another set of annihilation 
and creation operators. It is straightforward to see that to the linear order in h and modulo 
corrections proportional to didjTj^ipi), the solution with this initial condition is 



5x^^(k,x,r/) = _e^xx-ifc,.(x)-^kv(,.-.) . ^-^[^^x) - r,] H^%[q{v*{^) - v)] , (20) 

where q — k + 'kv. This is basically the Lorentz boost of the solution that one would find 

for r/^, = const. 

At smaU 77*(x) — rj, one has 6x2 oc [r^*(x) — r]]~^, i.e., the same behaviour as in (17). So, 
the phase perturbation freezes out: 



+ h.c. 



^ Rex(x,?7) J v^47r3/2(A; + kv) [ \ k ) 

(21) 

where wc again omit an irrelevant constant phase factor. Note that for 7]^ constant in space 
(and hence v = 0), i.e., to the leading order in /i, the phase perturbations are Gaussian 
random field with fiat power spectrum (4). The interaction with the radial perturbations 
makes the situation less trivial. 

The expression (21) serves as the initial condition for the evolution of the phase pertur- 
bations at the subsequent, nearly Minkowskian stage. We indicated in (21) that there is 
a correction of order dAjV*/^ (the factor is clear on dimensional grounds). The latter 
correction has been calculated in Ref. [11]; it will be irrelevant in what follows. 

2.4 Phase perturbations: order v"^ 

To calculate the tilt in Section 6, we will need the expression for 59 vahd to order v'^, but 
still to the first order in the gradient expansion of 577*(x) (i.e., corrections of order dAjV*/^ 
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arc still neglected). To this end, one observes [11] that to this order, the function (17) is no 
longer a solution to the field equation. One has instead 



Rex 



1 



jh[r)^{x) -r}] ' 



where 7 = {1 — v'^) ^1"^. Again using the analogy with the Lorentz boost, one obtains, instead 
of (20), 

(^X^^(k,x,,7) = ^n\\^{^,\^vn)M^-^--r<,^nm . ^|^g[^,(x) - ,7] Efl^Yiq{^,{^) - 77)] , 

where the indices || and T refer to components parallel and normal to v, respectively, the 
boosted momenta are 

?ll = 7(^11 + ^^^) , q^ = k^, g = 7(^ + ^11^^) , 

and, consistently neglecting the second derivatives of 5?7*(x), we have used 77*(x) = ?7*(0)— vx. 
In the limit q{ri^{x.) — r^) — )> one obtains the late-time expression for the phase, which can 
be written in a form, surprisingly similar to (21), namely 



h 



^47r3/27(A; + kv) 



gikx-ifcrj* (x) ^ 



1 + 



+ h.c. 



(22) 



The only difference with (21) is the factor 7 ^ = (1 — w^)^/^ in the integrand. 



3 Evolution at intermediate stage: order v 

As outlined in Section 1, our scenario involves the evolution of the phase perturbations S9 
from the hypersurface rj — 77* (x) to the hypersurface rj — rii — const. At this intermediate 
stage, the radial field stays at the minimum of the scalar potential, while the phase field 
is minimally coupled to gravity, and evolves in the sub-horizon regime. At time rii, the 
phase perturbations become super-horizon and freeze out again. The evolution of the phase 
must be nearly Minkowskian at this stage, otherwise its power spectrum would be grossly 
modified, see also Appendix A. So, the quantity of interest is 56'(x, r/i), and it has to be 
evaluated by solving the Minkowskian equation 

nse = {59)" - dAse = o . (23) 

The initial condition 5^(x, 77*(x)) at the hypersurface rj = ?7*(x) is determined by the dynam- 
ics at the conformal rolling stage. In this Section we perform the calculation to the linear 
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order in f, so the explicit expression is given by (21). The second initial condition is that 
the perturbation 59 is frozen out by the end of the conformal rolling stage, so that 



dN^e = at = r7*(x) , (24) 

where On denotes the normal derivative to the hypersurface 77 = 77* (x). As pointed out in 
Section 1, the case of interest is k{rji — rj^) 3> 1, so the evolution is long. 

3.1 Warm up 

It is instructive to begin with the unrealistic case 

77* (x) = 7740) - vx 

with constant v. This means that the Cauchy hypersurface is flat and boosted with respect 
to the cosmic frame. Let us consider the solution to Eq. (23) obeying the initial condition 
(cf. (21)) 

5^^k(x, 77*(X)) = e^kx-ifcr,.(x) ^ Q^^Q^ ^ . 

By going to the boosted reference frame back and forth, one finds that the solution, to the 
first order in v (and hence in /i), can be written as follows, 

5^k(x, 77) = Qi{^+kw)[^+^{r^-v*{m-ikvm cos[(A; + kv)(?7 + vx - 77*(0))] . 
Equivalently, 

The first lesson is that the solution is the sum of waves travehng along k and (almost) in the 
opposite direction; we will see in what follows that this situation is generic. Furthermore, 
for large enough {rj — rj*) the two terms in (25) have very different phases at given x, so 
their interference is negligible when integrated over k with any smooth function. The second 
lesson is that the wave moving along k has momentum k, while the momentum of the wave 
moving in the opposite direction is (k + 2A;v). We interpret this as the Doppler shift. Indeed, 
let us go to the reference frame (r, y) that moves with velocity v with respect to the cosmic 
frame, i.e., 

x = y-vT, 77 = T-vy 

(recall that we work to the first order in v). The Cauchy hypersurface rj = r)^{'x) corresponds 
to T = 7;*(0) = const, and the mode at this hypersurface is 



<^^k(y) 



_^ikx— ifc?7, (x) i(k+A;v)y _ —i(k+kv)r]» (0) 
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The last factor here is merely a constant phase, while the first factor describes the wave with 
momentum (k + kv) in the new reference frame. In the cosmic frame, this momentum gets 
shifted by —kv and A;v for waves moving along k and opposite to k, respectively. Hence the 
result (25). We will see that this situation is also generic: to the first non-trivial order in h, 
the main effect due to the intermediate stage is precisely the Doppler shift and the lack of 
interference between waves coming in the directions of k and — k. 



3.2 General formula and saddle point calculation 

The general solution to the Cauchy problem for Eq. (23) with the field and its normal 
derivative specified at hypersurface E is 



5e(x) = ^ dE'^ S^D^'^\x,y)d^5e(y) - 



dyi^ 



59{y) , (26) 



where D"^^* is the retarded Green's function of Eq. (23), x collectively denotes the coordinates 
(77, x), and the normal to the hypersurface is directed towards future. In our case the first 
term in the integrand is absent because of (24). We make use of the explicit expression 
(valid in the case x^ > y^ we are interested in) 

D^'\x,y) = ^S[ix-yy], (27) 

perform the integration over the radial variable and obtain for large (?7i— 77*) (see Appendix B 
for details) 

59{x) = f "^f" ^ rdrSe , (28) 



4:71 1 — nv 

where we still use the notation Vi = —dirj^,. Here n is unit radius- vector, integration runs 
over the unit sphere parametrized by n, and r = r(n) is the spatial distance that light travels 
from the hypersurface rj — 77* (y) to the point x — {r]i,x). It obeys the following equation: 

r ^r]i-r)^{x + nr) . (29) 

The function S9 — 66{r, n) in the right hand side of (28) is the field value at the Cauchy 
hypersurface, 

with 

y = x -|- nr . 

The formula (28) is exact for large r (for arbitrary r and general Cauchy data with non- 
vanishing dNS9, its generalization is Eq. (59) in Appendix B). 
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We now make use of (21) and obtain 

h r £k 

where / is the integral over unit sphere, 



^^(X' ^i) = / ^e^'^'^k • / + h.c. , (30) 



47r (l-nv)(A; + kv) 

with 

■0 — knr — A;77*(x + nr) = knr^i — (kn + A;)77*(x + nr) . (32) 

All quantities in the integrand of (31) (including v) are to be evaluated at y = x + nr. 
Corrections to the integrand are of order v'^ and dv/k. 

The exponential factor e*'^ in (31) is, generally speaking, a rapidly oscillating function 
of n, since ip is proportional to the large parameter kr. Therefore, the integral (31) can 
be calculated by the saddle point method, adapted to our problem. When performing the 
calculation, we have to keep in mind one point. Namely, even though we deal with soft 
modes in (5?7*(x) (with momenta p <^ k), the term kri^{'x + nr) in ip also gives rise to a 
rapidly oscillating factor, since r is large. So, we cannot neglect the second derivatives 9^77* 
in the exponent xjj. 

The saddle points are extrema of ilj{n), where n is a unit vector. To find them, let us 
formally consider n as an arbitrary vector, and 'ip formally as a function of this vector. Then 
the extremum on unit sphere is the point where dip/dn is parallel to n, i.e., 

^ = Xkrn (33) 
an 

with yet to be determined A (the factor kr on the right hand side is introduced for further 
convenience; in fact, Xkr is nothing but the Lagrange multipher). We use Eq. (29) to find, 
to the first order in v, 

dr 
on 

and, therefore, 

^ = [k+(kn + A;)v]r. (34) 
an 

We see that there are two saddle points, one near the unit vector k — k/k directed along the 
momentum, and another near (— k). These saddle points correspond to waves moving from 
the Cauchy hypersurface in directions opposite to k and along k, respectively, in accord with 
the discussion in Section 3.1. 
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The contributions of the two saddle points to the integral (31) are calculated in Ap- 
pendix C to the first order in v and dv. They sum up to 

^ = k b'""^ ~ ^""^^^^ + ''^^'^ ~ 4%)a,i;j+'^)] + e'^- (l - kv(-^)) } , (35) 

where 

-0+ = ip+i'x, k) = kr]i — 2A;77*(x + kr) , 

and superscripts (+k) and (— k) indicate that the corresponding quantities are to be evalu- 
ated at 

y(+) = X + kr (36a) 

and 

y(-) = X - kr , (36b) 

respectively. The terms in (35) marked by -|- and — come from the saddle points n fa k and 
n — k, respectively; they are analogs of the two terms in (25) (the factor {k + kv)"^ = 
k~^{l — kv) in the integrand in (21) was ignored in Section 3.1). Note that there is no 
symmetry between the two contributions; technically, this is because the dependence on Sr)^ 
is absent in the phase (32) for n = — k, but present for n = k. Note also that the saddle 
point value ip^ depends on x already to the linear order in h, while the second saddle point 
value ip- does not. This is precisely what we observed in Section 3.1: the momentum of 
perturbation corresponding to the first contribution in (35) is + dijj^/d:x. = k-|-2A;v, like in 
the first term in (25). Note finally that since we consider the case kr » 1, it is legitimate to 
neglect the correction of order d'^rf^/k — dv /k, indicated in (21), while keeping the correction 
of order rdv in (35). 

One more remark is in order. Our notation v*^^**^-* suggests that these quantities are 
functions of the direction of momentum only, i.e., that they are independent of the length 
of the vector k. This is true, but within our approximation only. The reason is that the 
horizon exit time rji is different for different k, so the arguments y^^^^ of v^^*^) depend on 
k through r = rji — rj^. This is irrelevant for us, since \i]i{k) — f]i{k')\ is at most of order 
1/k or 1/k' {in fact, it is even smaller, cf. Appendix A), so the effect we discuss is of order 
dv/k. Also, one may worry that the phases ip± depend on k through rji. This is irrelevant 
as well, for the following reason. When calculating the correlation functions of the field 59, 
one neglects the interference between the contributions due to the first and second saddle 
points, since the interference term oscillates in k as e^^^^ and is negligible when integrated 
with any smooth function of k. Then the factor, say, e*^''^ is merely a phase factor that can 
be absorbed into the redefinition of A^. In other words, x-independent phases cancel out in 
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the correlation functions of 66, so the dependence on k through rji does not appear. These 
observations apply to all calculations in this paper, so we neglect the dependence of rji on k 
in what follows. 

We conclude this Section by the discussion of the range of validity of our saddle point 
calculation. It follows from (32) that the relevant region of angular integration in (31) near 
each of the saddle points is A'j? ~ (kr)"^^'^. The saddle-point calculation makes sense if 
r7*(x + nr) does not change dramatically at this angular scale. Hence, by the saddle point 
method we can only treat the interaction of the phase perturbations with the modes of Srj^: 
whose momentum p obeys prA-d < 1, i.e., 

P^\ll- (37) 

The momenta p relevant for the statistical anisotropy do obey this inequality, see Section 4, 
while the requirement (37) restricts the angular scales at which we can reliably study non- 
Gaussianity. The latter point is further discussed in the end of Section 5. 



4 Statistical anisotropy 

We see from Eqs. (30) and (35) that the resulting phase perturbation (5^(x, r^i) is a combi- 
nation of two random fields, one associated with operators and AI. and another being 
5?7*(x). Let us discuss the two-point product 66{x)66(x.') averaged over the realizations of 
y4k and AI. for one realization of 5r;*, still to the linear order in h (in this Section we consider 
solely the resulting perturbations S9{'x,r]i) and omit the argument rji in the notation). As 
discussed in the end of Section 3.2, we neglect interference between terms with e^'^+ and e^'^- . 
Then the two-point function reads 



+ i/ ^e"""-"'-(l-2tv'-'))}, (38) 

where we made use of the fact that, to the first order in v, 

V'+(x, k) - V'+(x', k) = 2kv^+^\x - x') . 

Since we consider the long-ranged component of the field v, i.e., p k, we neglect the terms 
of order |x — x'| • dv. In particular, we do not distinguish between v(x' -|- kr) and v(x + kr) 
in the right hand side of (38). 

We now see explicitly that the actual momentum corresponding to the first term in (38) 
equals k -|- 2A;v, whereas the momentum in the second integrand equals k. To obtain the 
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standard form of the Fourier expansion, we change the variable to k = k + 2/cv in the first 
integral. To the first orger in h, the Jacobian of this change of variables is 




1 - 2kv(+'') - 2fc^5^ = 1 - 2kv(+'^) - 2divi^''^ ■ r(Sii - kh) , 
dki 



where we recalled that v'^^'^-' = v(x + kr). It is worth noting that the last term here cancels 
out the last term in square brackets in (38). So, omitting tilde over k, we obtain that for 
given realization of v(x), the power spectrum, with the correction of the first order in h, has 
the following form: 

Vse{k) = Po [l + k (^i^''^ - ^1"''^)] =Vo[l + g(k)] , (39) 

where 

Va = 



7r2 



is the power spectrum to the leading order in h (it is twice smaller than the power spectrum 
at conformal rolling stage after freeze-out of the phase perturbations; this is because the 
contributions of the two saddle points do not sum up coherently at rj = rji). Note that the 
non-trivial term in (39) depends on the direction of momentum k. Note also that the power 
spectrum (39) is symmetric under k — > — k, so the two-point function (38) is invariant under 
X ■<-> x', as it should. Low angular harmonics of v^^'*\ viewed as a function on unit sphere 
in momcntTim space, take certain values in our patch of the Universe. Hence, they induce 
statistical anisotropy; in particular, the lowest multipole of the expression in the right hand 
side of (39) (quadrupole) gives rise to the power spectrum of the form (9), (10). 
In more detail, the right hand side of (39) contains all even multipoles, 

g(k) = J]?i^i^^(k), (40) 

Im 

where Yim are spherical harmonics. Making use of the definition v^='='*) = v(y (=*=)), where y^='=) 
are given in (36), we find for Z 7^ that the multipole coefficients are given by 

where we omitted an irrelevant k-indcpcndcnt phase. It is worth noting that for low multi- 
poles, the relevant range of integration over p is roughly p ~ r^^: at larger p the integrand 
rapidly oscillates, while at smaller p the expression in the inner integrand in (41) decays as 
while according to (18) the amphtude of Srj^{p) behaves as ^/Ps^ oc p~^. At large I, the 
relevant momenta are of order p ~ Ir"^. Thus, our approximation p <^ {k/rY^^ is justified 
at least for low multipoles. 
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The calculation of the variance of qim is performed in much the same way as the calculation 
of the CMB anisotropy multipoles, see, e.g., Ref. [28]. This is done in Appendix D with the 
result 

Note that we use different normalization here and in (10). To establish the correspondence, 
we calculate the angular integral of the variance of the quadrupole term in (10): 



while the same integral of the quadrupole term in (40) is given by 

/2 2 
m=-2 ' 



m=-2 



Hence the extra factor 75/87r in (11) as compared to (42). 



5 Non-Gaussianity 

The statistical anisotropy is an appropriate notion for describing the effect due to the varia- 
tion of v^^^^ over large angular scales in momentum space. On the other hand, the effect of 
fluctuations of v^^**^) at small angular scales is naturally interpreted, we believe, in terms of 
non-Gaussianity. Indeed, in the latter case it makes sense to treat v'^^'^) as genuine random 
field and perform averaging over its realizations, having in mind multiplicity of patches in 
the k-sphere. 

It is worth noting that even though we are going to consider v^^"*^ at small angular scales 
Ai? in momentum space, the relevant momenta p of the field ^Ty* are still small, p ~ (rA'j?)"^. 
So, our approximation p <^ (k/rY^^ is still valid, provided that Ai? is not very small, see the 
discussion in the end of this Section. 

Let us consider higher order correlation functions of 66 (x.) (we again omit the argument 
f]i in this Section). Since this field has the general structure (30), where / does not contain 
the operators Ak, AI_, the three-point function vanishes identically. For calculating the non- 
Gaussian part of the four-point function, the expression (35), valid to the first order in v, is 
sufficient. Proceeding in the same way as in the beginning of Section 4, we obtain 

{59{k)59{k)69{k')S9{k')) = -l^-L^p^ s{k + k)S{k' + k') 

X 1 + GNcik k') + GNG{-k, k') + GNcik -k') + GNci-k, -k') 

+ (k ^ k') + (k ^ k') , (43) 
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where the non-Gaussianity is encoded in 

G.aik kO = {k {v^' - vr^) ■ k[ [v\^^'^ - .;-^'))) . (44) 

Fluctuations of v^^''^ at small angular scales in momentum space show up when k and k' 
are either nearly parallel, or nearly antiparallel, the latter case being related to the former 
by the interchange k •H- k. So, it suffices to consider nearly parallel k and k', i.e., 

|k-k'| < 1 . 

Since the power spectrum of the random field v(x) is flat, the leading term is logarithmic in 
|k-k'|. 

The expression in (44) involves the combination 



.(+k) 



Therefore, the integral over momenta of the field v is cut off in the infrared at p ~ r^^. 
We cannot quantitatively treat the modes of these momenta anyway, since they are plagued 
by cosmic variance. So, we consider modes with p > r~^, recall the expression (18) for the 
power spectrum of ^Ty* and write 



The angular integral here is straightforwardly evaluated. We make use of the fact that 
k(k - k') = 0((k - i!f ) and obtain 

^ 9/i^ C dx f sin a; 

47r2 X>|k-ic'| x'^\x 

where x = rp|k — k'|. This is a logarithmic integral, and in the leading logarithmic approxi- 
mation we immediately get 

3h^ const 

Gng = -r^ log ^ . 

47r^ |k-k'| 

The constant here is of order 1; it cannot be rehably calculated, since the contribution of 
the region p ~ is undetermined because of the cosmic variance. Finally, we notice that 
the right hand side of (44) is symmetric under k ^ — k, so the four terms in (43) give equal 
contributions. Thus, the four-point function at |k — k'| <C 1 is 



{59{^.)5e{i.)69{^.')59{i.')) = ^,^vl 5(k + k)5(k' + k' 



l + FjvG(k-k') 
+ (k ^ k') + (k ^ k') , 
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where 

7-. /f f'\ const 

FivG(k - k') = ^ log ^-3^ ■ (45) 

We conclude this Section by noting that our analysis is valid provided that we can treat 
the range x = rp|k — k'| ~ 1 within our approximation p <^ {k/r)~^^^, see (37). So, the 
logarithmic behaviour (45) persists until |k — k'| > (A;r)~^/^. At even smaller angles between 
k and k', the function F/vG'(k — k') most likely flattens out. The logarithmic behaviour is 
definitely absent for |k — k'| < (kr)^^, since momenta p higher than k do not contribute to 
the effect. These observations suggest that the value of r is potentially measurable. 

6 Tilt 

Once the interactions of the phase field with the radial one are not neglected, the power 
spectrum of the phase perturbations obtains a small tilt. The reason is that for larger k, 
there are more modes of 5?7* with p < k which affect the properties of the phase perturbations. 
We will see that the effect is logarithmic because of the flat spectrum of v. 

To this end, let us come back to the two-point correlation function {S9{yi)S9{x.')). Even 
though the integral (28) is again saturated near n = ±k, the saddle point calculation like that 
performed in Section 3.2 is no longer appropriate, since we are going to consider all modes of 
5?7* of momenta p <^ k and not necessarily very large wavelength modes obeying (37). The 
problem is not notoriously difficult, nevertheless, since we are interested in logarithmically 
enhanced effect. Imagine that one calculates {59{^)59{^')) by expanding in (5?7* the complete 
expression (28), with 59 in the integrand given by (22). In principle, large logarithms could 
come from the expectation values {Sr]^ ■ didjSr)^) and {vi ■ Vj). We reiterate, however, that the 
overall time shift is irrelevant for our problem, so the terms of the former type do not appear 
explicitly (for the same reason, there are no terms involving correlation functions of Sr)^ with 
itself and with v, which would yield power law corrections). Thus, it is legitimate to ignore 
the correction of order didjdr],, in (22). Moreover, we can formally consider the velocity v 
in (22) as a constant which is independent of spatial coordinates. So, we effectively deal 
with the Lorentz-boosted hypersurface = rjtf{y±) — v(y — y±), where y± = x ± kr. The 
qualification here is that the velocity is to be evaluated at y = y±, and that v(y-i-) is a 
non-linear function of 5r]^, since, according to (29), y± depend on 5r]^ through r. Another 
qualification is that when calculating the power spectrum Vse at momentum A;, we have to 
impose a restriction p < A; on the momentum p of modes of the field ^77*. 

Since we treat the velocity v as constant in space, we can obtain the solution to the 
Cauchy problem explicitly, in a way similar to that of Section 3.1. However, we need the 
solution to the second order in v. The initial condition for the Minkowskian evolution is thus 
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given by (22). Let us define the Lorentz-boosted coordinates z and r: 



z\\ ^ -f{x\\ + vrj) , zt = xt, t ^ -f{rj + vx\\) , 

where || and T refer to components parallel and normal to velocity. Then the initial data 
are specified at the hypersurface t = t± = const, and 8-^69 = at this hypersurface. We 
re-express x and rj in terms of z and r and insert them into (22). Omitting the overall 
phase factor independent of z that cancels out in the two-point function, we write the initial 
conditions as 

/I d^k 
— — — e'^iMk^ + h.c. , dje = , (46) 
7(/c + kv) V2A; 

where 

QW = lih + kv) , qT = kr . 
The solution to the massless field equation in Minkowski space with this initial condition is 

/I d^k 
-j^e- cos + /,.., (47) 

where q — ^{k -\- k\\v). This solution again describes two waves propagating in opposite 
directions, which do not interfere at 77 = 771. Let us consider the two waves separately. 
At time 771, we have for the first wave, moving in direction opposite to k, 

where </? is a phase, irrelevant for the two-point function of 59. So, the actual momentum is 

k\\ = 7^(A;|| + 2kv + k\iv^) , k^ = k^ . 

Note that to order v'^ we have A;|| = ^2v{k\\ + '^kv), where 72^ — — ^v^)"^^"^ is the Lorentz- 
factor for velocity 2v. Hence, k again differs from momentum k by the Lorentz-boost with 
velocity 2v. We recall that Ax^d^k/ \/2k is Lorcntz-invariant, and obtain for the contribution 
of the first wave at point x (again omitting the phase factor, irrelevant for the two-point 
function) 

(5^(x)oc / ^_^e^kx^^4t + /i.c., 
J 7(/c + kv) 

where k\\ = 'y2v{k\\ — 2kv), k = '~f2v{k — 2k\\v). Expanding in v to the second order, we obtain 
the following form of the first contribution to the power spectrum 
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Here the term in parentheses comes from {69^^^ 59^^^), where 59^^^ is the zeroth order phase 
perturbation and 69^^^ the correction (that includes hnear and quadratic terms in v), while 
the last term in the right hand side is due to the correlator {69^^^S9^^^). We see that the 
explicitly quadratic terms cancel out and find (at this point we can set k = k) 

Vse<xl + 2ki{vi) . 

We now recall that the velocity is to be evaluated at the point y+ = x + kr, where r = 
Vi ~ ^*°^ — 577*(x + kr), so that 

Vi{y+) ^Vi[x + k{rji - rji^'^)] - djVi ■ kjSrj* . 

The expectation value of the first term on the right hand side vanishes, while the second 
term gives 

{vi{y+)) = -kj J -^PiPj'^Sv* (P) = -^i-^ \og{kr) , 

where we recalled that the relevant range of momenta is <^ p <^ k. Thus, the contribution 
due to the first wave has the form 

3/^2 

Vse(xl-—log{kr) . (49) 
Let us now consider the second wave that moves along k. We have at time rji 

SO the actual momentum is equal to k. Hence, the contribution of this wave is 

S9{^) oc / -^—^e^^A^^ + h.c. . 
J 7(/c + kv) ^/2k 

Proceeding as before, we obtain the contribution of this wave to the power spectrum, 

Vseocl- 2ki{vi) , 

where the velocity is to be evaluated at the point y_ = x — kr with r = rji — 77^°^ — 5?7*(x — kr). 
The resulting contribution again has the form (49), so we conclude that the shape of the 
entire power spectrum is given by (49). 

The result (49) shows that the power spectrum of 59, and hence of the adiabatic pertur- 
bations, is tilted. If this is the only reason for the tilt, the scalar spectral index in our model 
is equal to = 1 — As pointed out in Section 1, however, there may be other sources 
for the tilt, so we cannot insist on attributing the entire scalar tilt to the effect discussed in 
this Section. 
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To end up this Section, we sketch an alternative way of calculating the correction to the 
power spectrum Vse- One makes use of the exact formula (28) with 59 in the integrand 
given by (22) and evaluated at y = x + nr, where r obeys Eq. (29). The dependence of the 
integrand in (28) on the integration variable n is fairly non-trivial, since the vector n enters 
the argument of 77* both explicitly and through r(n). We know, however, that the integral 
(28) is saturated in regions near the two points on unit sphere, n = k and n = — k. Consider 
the first region for definiteness. The idea is to write 

?7*[x + nr(n)] = r7*[x + kr(k)] + j^T^Jx + nr(n)] - (5r7*[x + kr(k)]| , 

r(n) = r(k) - |(5r7*[x + nr(n)] - (^?7*[x + kr(k)] | , 

express these functions iteratively through 

(5?7*[x + nr(k)] - 5r]^[y. + kr(k)] 

and systematically expand the integrand in (28) in a series in the latter quantity, up to 
quadratic order. Then one has to deal with angular integrals, in which the integration 
variable n enters either in combination e*^"''^'*) or via e*^"''^'*^5r7*[x-|-nr(k)] (the integral with 
8ril is trivial after ensemble averaging). The former integral is straightforwardly evaluated 
by the saddle point method. To evaluate the latter integral, one writes 5?7*[x-|- nr(k)] in the 
Fourier representation and arrives at the angular integral with e**^''"''^)"'"^''^ where p is still 
the momentum of a mode of 57]^. The latter integral is again evaluated by the saddle point 
method; the rest of the calculation is straightforward. 

We have performed the calculation of the power spectrum in this way; it is tedious, but 
does yield the result (49). 
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Appendix A. Scalar field perturbations in a contracting 
Universe with super-stiff matter 

In this Appendix we discuss the free propagation of masslcss scalar field $ minimally coupled 
to gravity in contracting Universe filled with matter whose equation of state is super-stiff, 
p — wp, w^l. Since the phase 9 behaves precisely in this way after freeze out of the radial 
field 1 01, our discussion applies directly to the situation studied in this paper. Our point is 
to show that in the limit w — > oo the propagation is effectively Minkowskian all the way 
down to a{r]) — >■ 0. 

For constant w, the scale factor evolves in conformal time as follows, 

a^\r]f , 77<0, 

where 

^ l + 3w 

In terms of the field a — a^, the field equation reads 

a" + - -a = a" + /tV + ^ii^a = . (51) 

For large w and hence small ^, the last term in the left hand side of Eq. (51) is negfigible 
before the horizon exit time, 1]^^ ~ —^/P/k, while there is simply no time to evolve even 
in Minkowski space in the time interval (r^exjO). This is why one can make use of the 
Minkowskian evolution to evaluate the value of the field $ as 7y ^ 0, i.e., deep in the 
super-horizon regime. 

To substantiate this claim, let us consider the Cauchy problem similar to that discussed 
in the main text. Namely, let the initial value $i be specified at 77 = 77* = const with 
|?7*| 3> k~^, and another initial condition is = at 77 = 77*. Let us compare the values of $ 
obtained at 77 = by solving the Minkowskian evolution equation = and by evolving 
the field according to Eq. (51). The Minkowski evolution gives ^^'^'^^{jj) — ^iCOsk{r] — r]*), 
so that 

$^^"^(77^0) =$iCosA;77* . 
The solution to Eq. (51) with the above initial conditions imposed at |77*| » is 

(7(77) = Hr]*fJ^\ [uHi'\-kr]) + u*Hi'\-kr])] , 
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where v = 1/2 — /3, 

^ ^ lgife»?*+if (1-/3) 

2 

and ifi^^'^^ are Hankel functions. The asymptotics of $ = cr/a as 77 for ^ < 1/2 is 
^[r) ^ 0) = $i cos /cr^* ' ' ' 



2 J \ 2 J r(l/2) 

We see that the Minkowski an result indeed coincides with the exact one in the hmit w — )■ 00, 
i.e., /3 — 7> 0. The main effect for finite but large w is the induced tilt in the power spectrum. 
The phase 7r/3/2 is irrelevant, as it cancels out in the correlation functions. 



Appendix B. Derivation of the formula (28) 

In this Appendix we consider the Cauchy problem for Eq. (23) with initial data specified at 
the hypersurface 

f(y)-v-v*{y)-o, (52) 

where y denotes a point with coordinates = (77, y). We simplify the notation and use 9{x) 
instead of S9{x). 

Let 9{x) be the solution to the D'Alembert equation (23), such that 9{y) and dN9{y) 
coincide with the Cauchy data 9{y) and dN9{y) at the Cauchy hypersurface (hereafter 
denotes the normal derivative). Let us introduce 

9{x) = 9{x) ■ Q[f{x)] , 

where © is a step function. Then 

D9 = d,9d^f ■ S{f) + d,[9d^f ■ 5(f)] 

and, therefore, 

e{x) = J d'y i^D"'\x,y) d,e{y) d^f{y) ■ 5[f{y)] - 

(53) 

where we omitted tilde over 9 in the right hand side, since the integration runs over the 
Cauchy hypersurface. The second term in the integrand is obtained by integration by parts. 
The formula (53) is nothing but the general formula (26), and 8^98'^/ oc 8^9. 

In the case of interest, the normal derivative vanishes at the Cauchy hypersurface, and 
the first term in the integrand in (53) is absent. We make use of (27) and write 

A_Dret^^,y)^-hx,-y,)S'[{x-yr]. 



d 

dyi^ 



D^^\x,y) 



Oiy)d''f{y)-S[f{y)] \ , 
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We use the explicit form (52) of f{y), integrate over rj in (53) and obtain for x = (r7i,x) 

9{x) = i y d'y [771 - 77,(y) + v(x - y)] ^(y) 5' ([r^i - v^y)]' - (x - y)^) , (54) 

where Vi — —dir]^{y) and ^(y) = ^[y, 77*(y)] is the field value at the Cauchy hypersurface. 
We now introduce the integration variable r via y = x + r, write r = nr, where n is a unit 
vector, and cast the integral (54) into the following form: 



e{x) 



— J (in„r^(ir [771 — 77*(x + nr) — nvr] ^(x + nr) 5' ([771 — 7y*(x+ nr)]^ — r^) . (55) 



Here v = v(x + nr). Finally, we make use of the identity 

1 d 

{[Vi - ^*(x + nr)]^ - r^) = -— — ■ {[rn - r?*(x + nr)]^ - r^) , 

^ \ /J / 2 |r — nv[?7i — 77*(x + nr)J} ar ^ v /j / 

which is obtained by evaluating the derivative over r of 5 ([771 — 77* (x + nr)]^ — r^). Since 
r 7^ at the Cauchy hypersurface, we can integrate over r in (55) by parts. We also use the 
fact that 

5 ([771 - 77.(x + nr)]2 - r^) = —^——^5[r - r(n)] , (56) 
where r(n) is the solution to Eq. (29). We get 



e{x) 



1 f d f r'^ \ 

- / dilndr — — — ■ -r-[77i - 77*(x + nr) - nvr] 9{x + nr) 

■K J 9r \2 {r - nv[77i - 77*(x + nr)J j / 



X — -S\r - r(n)l 

2r(l - nv) ^ ^ 



The integration over r is now straightforward, and we obtain after some algebra (note the 
cancellation of the terms with derivative 5v(x + nr)/dr) 



e{x) = ^J da 



9 + rdr9 

1 — nv 



(57) 



where in the right hand side one has 9 = 6'(y,77*(y)) with y = x + nr. Let us emphasize 
that (57) is the exact result for the Cauchy problem with 8^9 — 0. At large r, the second 
term in the integrand dominates, and we arrive at the formula (28) used in the text. 

For completeness, let us derive the general formula for the solution to the Cauchy problem 
with non- vanishing 8^9. With the Cauchy hypersurface defined by Eq. (52), the derivative 
along the unit normal is given by 

8N9 = -f8^9 8>'f , (58) 
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where 7 = (1 — v^)"^/^. This expression can be obtained by performing local boost 

dr — ^{dr] + vdx) , etc. 
Then r is the time coordinate along the normal, and 

dj,e = dj = ^{dr,e - vAO) , 

which is precisely (58). Making use of (58) and (27) we write the first term in (53) as follows, 

/ A^<^(k-^*(y)r-(x-y)^) l^dr^e. 

We proceed as before, again use (56) and obtain for this term 

If r 
Att J 7(1 — nv) 

Thus, the complete expression for the solution to the Cauchy problem is 

9(x) = — [ dn^ 9+ — - — r ( drO + \/l^dNe] . (59) 
Att J [ 1 - nv \ /J 

The notations here are the same as in (57). 

Appendix C. Details of saddle point calculation 

Saddle point n « k 

To find the saddle points of the integral (31), we solve Eq. (33) with dip/dn given by (34). 
To the linear order in h, the first saddle point is 

n+ = k + 2[v-k - (kv)] 

with 

A = 1 + 2 kv . (60) 

Let us evaluate the contribution to the integral (31) coming from the saddle point region 
near n_|_. Let (p be angular coordinates in the frame with the third axis along n_|_. Then 

n = n+ + n(^) + n^^^ , 

where n^^^ and n*^^) are of the first and second order in i?, respectively, 

n*-^-* = (sin'(9cos(/9, sin?? sin(^, 0) , 
n(2) = (0,0, cos?? - 1) . 
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We have 

^(n) = ^(n+) + ^/;(^) , 

where 

and the derivatives are evaluated at n = n+. The first derivative is given by Eqs. (33) and 
(60), while to the linear order in v and dv (i.e., linear order in /i), the second derivative is 



driidrij 



r{kiVj + kjVi) + r (kn+ + k)diVj . 



The angular integral is now straightforwardly evaluated (one first integrates over i? near 
— with weight '&d'&, then expands in v and dv and integrates over (p), and to the linear 
order in h one finds 

The pre-exponential factor in (31) is to be evaluated at n = n+. Collecting all factors, we 
get the contribution of the first saddle point (to the first order in h): 

^2 A; + kv 

Note a non-trivial cancellation between v-dependent terms in the pre-exponential factor. 
Finally, we recall that 

■0(n+) = kn+?7i - (kn+ + k)r)^{-x + n+r) = krji - 2kr]^{x + kr) , 

where we still work to the linear order in h. Since Srj^, and v are already of order h, their 
argument is merely y^"*'-' = x + kr. In this way we arrive at the first term in (35). 

Second saddle point 

The second saddle point is precisely at 

n_ = -k 

(this is exact result valid to all orders in v). At this saddle point we have 

■^(n.) = —kfji . 

The same calculation as above gives for the contribution of the second saddle point 



2 A; + kv 

So, the second term in (35) is obtained in a very straightforward way. 
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Appendix D. Multipoles of statistical anisotropy. 



The field 5rj^{x) is an isotropic Gaussian field. Therefore, the multipole coefficients in (40) 
are independent, 

illmQl'm') — Ql^ll'^mm' ■ 

We make use of the expression (41) and calculate the sum J^^d P). Since {Srj*{p)Sril{p')) oc 
S{p — p'), this sum has the following form: 

E(l^-I') = / d^^sY.l'iimip)]'- (61) 

m in 

The integrand here is independent of the direction of p and therefore can be calculated in 
any reference frame. To simplify formulas, we choose, somewhat loosely, a reference frame 
one step earlier, in the inner integral in (41), so we calculate qim{p) in a p-dependent frame. 
This procedure is legitimate as long as one calculates the sum in the right hand side of (61). 
We choose the spherical frame with p directed along the third axis and write 



1 

= -^(5„oV(2/ + l)7r J dtPiit) ■ pt ■ {e'P'' - e-'^') , (62) 
-1 

where Pi are the Legendre polynomials, i) is the angle between the momenta p and k and 
t = cost?. Since the integrand in (62) is symmetric under t ^ —t (this is a consequence 
of the symmetry of the power spectrum Vd0(k) under k — >■ — k, see (39)), odd multipoles 
vanish. In what follows we consider even 1^0. 

The standard way of calculating the integral (62) is to make use of the expansion of the 
oscillating exponent in Legendre polynomials, 

oo 



l'=0 



where ji are spherical Bessel functions. We make use of the normalization of the Legendre 
polynomials. 



^dtPi{t)Pi,it) ^ ^j^Sw , 



and recurrence relation 

l'Pi,^,{t) + {l' + l)Pi,+,{t) 



tPi'it) = 



21' + 1 



31 



Then the integral (62) is straightforwardly evaluated, 

qimip) = 2dmo^l^^ji^p[{l + ^)ji+i{y) - hi-1 
where 

y^rp. 

We now insert this result into (61), recall that the power spectrum of 61]^ is given by (18) 
and get 



i-i{y)\^ ■ (63) 



Finally, we recall the relationship between the spherical and conventional Bessel functions. 



and perform integration by using 



My)My)y-'dy^ r(A)r (^+-/+^) 



After straightforward algebra this yields 

v^,, ,2- 3h^ 21 + 1 , ^ 

V'"^'^^V (Z-l)(Z + 2) ' 

or, equivalently, the quoted result (42). 

It is worth noting that the relevant integration region in the integral (63) is y = pr ^ I 
(the spherical Bessel function ji{y) is exponentially small at y <^ Z and decays as y~^ at 
y ^ I). This means that our approximation p <^ {k/rY^^ is justified for /cr » 1, unless I is 
very large. 
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